Limits to multipartite entanglement generation with bosons and fermions 
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Many-photon interference in linear-optics setups can be exploited to generate and detect multi- 
partite entanglement. Without recurring to any inter-particle interaction, many entangled states 
have been created experimentally, and a panoply of theoretical schemes for the generation of various 
classes of entangled states is available. Here, we present a unifying framework that accommodates 
the present experiments and theoretical protocols for the creation of multiparticle entanglement via 
interference. A general representation of the states that can be created is provided for bosons and 
fermions, for any particle number, and for any dimensionality of the entangled degree of freedom. 
Using the framework, we derive an upper bound on the generalized Schmidt number of the states 
that can be generated, and we establish bounds on the dimensionality of the manifold of these states. 
We show that - at the expense of a smaller success probability - more states can be created with 
bosons than with fermions, and give an intuitive interpretation of the state representation and of 
the established bounds in terms of superimposed many-particle paths. 



I. INTRODUCTION 

The preparation of entangled states of photons [1] has 
been the subject of a great collective effort, which was 
driven by fundamental scientific questions such as the 
violation of local realism [2HI] as well as by the perspec- 
tive of technological applications, including, e.g., quan- 
tum cryptography [5J and quantum computing 6 . The 
control over the process of spontaneous parametric down- 
convertion (SPDC) started this successful journey in 
1988, when a two-photon entangled state was created 
[7J. This achievement was followed, in 1999, by the cre- 
ation of an entangled tripartite state [5]. Step-by-step- 
wise, the number of photons was further increased, such 
that entanglement in four-photon [9], five-photon [ID] , 
six-photon [TTMl4j . and, eventually, eight-photon states 
[151 116] was demonstrated. While many experiments can 
generate only one particular state, there are also setups 
that can tune the final state through several classes of 
qualitatively different types of multipartite entanglement 
[17] by varying parameters in the scattering setup [T5] . 

The creation of an entangled state of photons relies on 
the independent propagation of the particles in a linear- 
optics setup and on the subsequent selection of successful 
final states with one particle in each spatial output mode, 
i.e. on post-selection. No direct interaction between the 
photons takes place, in contrast to other possible car- 
riers of entanglement such as trapped ions [19] . It is, 
in contrast, the indistinguishability of the photons which 
provides the seed for quantum correlations. Therefore, 
the physical mechanism - the coherent superposition of 
many-particle paths - is rather independent of the prac- 
tical experimental implementation. For example, also 
the indistinguishability of cold atoms can be used for 
interaction- free entanglement generation [20] , 

Whereas in the mentioned examples for multiphoton 



entanglement, quantum correlations are encoded in the 
polarization of the photons, also other degrees of freedom 
can be used to carry entanglement, such as the time-bin 
[2TH23] , the path [7j [24] , or the orbital angular momen- 
tum [251428] . Since degrees of freedom other than po- 
larization are not restricted to two dimensions, their use 
renders the creation and detection of entanglement in 
higher dimensions possible. In particular, genuine qutrit- 
cntanglement was realized with energy-time entangled 
photons [53]) and the encoding of a 20-dimensional de- 
gree of freedom in the time-bin of was demonstrated with 
a single photon [3D] . 

Besides the variety of entangled states that has been 
created experimentally with photons by now [TJ 1131 1141 
ITS] [31] . manifold further strategies for the creation of 
large classes of entangled states were proposed [3"2H3"8"] ; 
but although all these schemes share the same physical 
ingredients, no common framework permits the direct 
comparison of these approaches. With the rapidly in- 
creasing complexity of many-particle interference setups 
with growing particle number, a general framework is 
highly desirable for the design of setups with high success 
rate and low number of components, and for the system- 
atic comparison of competing approaches to multipartite 
entanglement generation [37, 138] . A common language 
for a unified description of such existing schemes, and 
for possible new strategies is our first desideratum in this 
article. 

From the conceptual perspective, we also lack an un- 
derstanding for the intrinsic limitations of setups that 
rely on many-particle interference. At first sight, it seems 
that a very large set of entangled states of photons of 
nearly arbitrary particle number and dimension can be 
created in the laboratory. It was, however, also recog- 
nized that a toolbox composed of linear optics and post- 
selection alone is by far not universal, i.e. a setup that 
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relies uniquely on the indistinguishability of particles im- 
poses severe restrictions on the implementable quantum 
operations. Consequently, several no-go theorems were 
formulated for few-particle systems. For example, a Bell 
measurement that can discriminate all four Bell-states 
is impossible with linear optics |39j . Similarly, teleporta- 
tion with 100% success cannot be achieved using photons 
[40) . On the other hand, the manipulation and char- 
acterization of qudits, i.e. higher-dimensional degrees of 
freedom, is severely constrained: The Schmidt-rank of 
any measurement operator on two qubits is bound by N, 
when N — 2 auxiliary qudits are available [UJ . 

While the characterization and the classification of 
general multipartite entangled states is rapidly progress- 
ing [4"2"H4"T] . it becomes desirable to understand the gen- 
eral physical principles that limit the creation of mul- 
tipartite qudit-likc entangled states. Also for this en- 
terprise, a common framework that describes generic 
many-particle interference-based setups is a mandatory 
requirement. Such framework should, both, describe the 
physical processes that lead to a certain final state, and 
offer a representation of this state that permits an in- 
vestigation with the tools of quantum information the- 
ory, like entanglement measures. In other words, such 
framework should bridge the gap between the physical 
many-particle scattering process and the abstract quan- 
tum state of quantum information theory. 

Here, we set the stage for a systematic investigation of 
the aforementioned problems. We propose a treatment 
of general entanglement-generating setups which accom- 
modates virtually all possible schemes that use identical 
particles, where the particles carry (typically internal) 
degrees of freedom that may or may not be entangled. 
Thereby, all existing schemes for photons can be under- 
stood in equal terms: The created state emerges as a co- 
herent superposition of many-particle paths [48 , which 
can be traced back to the state representation. This state 
representation allows us to investigate systematically the 
capabilities of interference-based setups and to character- 
ize the set of states that can be created, with the tools 
of quantum information theory. In particular, we find a 
general combinatorial bound for the maximal generalized 
Schmidt number [3!5] for multipartite entangled states of 
qudits, which generalizes a previously known theorem for 
two particles [UJ . This bound constitutes an intuitive re- 
striction of the set of states that can be generated, and 
a physical characterization of the intrinsic limitations of 
interference-based entanglement generation. It applies 
equally for bosons and fermions, for which, however, the 
dimensionality of the set of states that can be created is 
different, as we will also show. As a consequence, more 
states can be generated with bosons than with fermions, 
at the expense of a lower success probability of entangle- 
ment generation for bosons. 

We first establish the general experimental capabilities 
that are available today, and those which may become 
available in the near future, in Section [TTJ On that ba- 
sis, we formulate the requirements that a versatile model 



should fulfill, and develop a framework that satisfies these 
demands. An explicit representation of the achievable 
multipartite states is then derived in Section |III| which 
allows the formulation of a general bound on the Schmidt 
number of the generated states. We then show in Section 
|IV| that this bound is not the only restriction that applies 
to interference-based setups, since the number of dimen- 
sions of the manifold of states that can be created with 
bosons and fermions establishes another limitation, with 
a tighter bound for fermions than for bosons. Concluding 
discussions are given in Section [V] For illustration, and 
in order to familiarize the reader with the methods and 
with our framework, we schematically reproduce several 
state-of-the-art experiments and theoretical schemes in 
Appendix [A] The derivation of the dimensional bounds 
is presented in Appendix [B] 

II. UNIFIED FRAMEWORK FOR 
ENTANGLEMENT GENERATION 

A. Physical situation and versatility requirements 

The examples for entanglement generation that were 
mentioned in the introduction are all based on the same 
principle mode of operation: Several photons prepared 
in a given initial state propagate through a linear-optics 
setup. Post-selection is performed on events at the out- 
put modes, and only events with exactly one particle in 
each spatial mode arc taken into account. Given such a 
post-selected, successful event, the internal states of the 
particles have then the desired entanglement properties. 
While photons are destroyed in the moment they are de- 
tected in an optical mode, the presence of massive parti- 
cles in a target mode can be ensured by some measure- 
ment that leaves the particle and the entangled degree of 
freedom it carries unaffected, and the entangled particles 
may be used for further processing. Besides entangle- 
ment generation via propagation and post-selection, also 
the projection onto entangled states and the detection 
of entangled states is possible with techniques that are 
based on the same principles [UH [5U] . 

Given the rather broad variety of existing schemes and 
successful experiments, our dedicated attempt to formal- 
ize the state-of-the-art under one housetop should fulfill 
the following requirements: 

• Any number of particles N is allowed. 

• Entanglement is established between arbitrary d- 
dimensional degrees of freedom. 

• The coherent manipulation of the entangled degrees 
of freedom is possible without restriction. 

• The initial state may or may not be entangled. 

• Bosons as well as fermions may carry the entangle- 
ment . 



3 



• The creation of entanglement as well as the pro- 
jection onto entangled states are contained in the 
description. 

In contrast to the polarization of photons, which can 
be manipulated essentially without restrictions, it is not 
immediate that all coherent single-particle operations or 
state-dependent beamsplitters can be implemented for 
any other degree of freedom [5T] . Depending on the cho- 
sen physical carrier of entanglement - e.g. polarization, 
time-bin, path, or orbital angular momentum -, it may be 
possible that, e.g., operations that discriminate the inter- 
nal degrees of freedom are difficult to implement [HTJ 155] . 
It is therefore important to pay particular attention to 
the modeling of different types of setups with distinct 
capabilities, as we will do in Section |IIE| 



B. Many-particle states and their interpretation 

A central building block of our present approach is 
the description of the physical many-particle scattering 
process and its relationship to the entangled states it 
generates. We emphasize the conceptual difference be- 
tween the abstract, quantum-information interpretation 
of a state (the "iV-qudit-state" ) , and its physical realiza- 
tion, which we will need to work with in parallel. 

For the understanding of an experiment, it is useful 
to describe physical states of photons or other identi- 
cal particles in second quantization; we will denote these 
states by \^> 2q ). In this representation, we can formulate 
the physical many-particle scattering process, whereas 
the definition of an entanglement content requires the 
identification of a subsystem structure which is most ad- 
equately done in first quantization |53) . Within a clear 
tensor structure that is imposed by physical observables, 
the usual means to characterize entanglement, such as 
monotones and measures |54j . can then be applied on 
the first-quantized state 

Since the dynamics is best described in second quan- 
tization, whereas first quantization is more appropriate 
for the description of entanglement properties, we will 
specify a map, Q a , between both representations such 
that I* 29 ) = fl a |$ l9 ) . We emphasize that the conceptual 
problem of the assignment of entanglement to a state of 
identical particles 55-57. is circumvented by the present 
treatment, since states of identical particles \^ 2q ) are 
identified with states of distinguishable subsystem com- 
ponents Before we can define this map, let us 
briefly review the essential properties of quantum states. 



1. State representation 
An JV-qudit-state, 

d d 

l $l9 ( c )> = £- E c n _ 3N \ n ,..., ]N ), (l) 

jl = l JJV = 1 



is fully specified by the d N components of the coefficient 
tensor 



c e <C d ® <C d ® • • • ® <D d . 



(2) 



The subsystem structure in Q and ^ reflects the 
paradigmatic situation of N parties that each control a 
(i-dimensional subspace. Therefore, the usual notions of 
separability and entanglement can be applied on |$ l9 (c)) 
Q21[SB1[SH]. For instance, the state |* 1<J (c)) is fully sep- 
arable (or iV-separable) if Cj u . 
product, 



can be written as a 



„« J 2 ) 



(3) 



The iV-qudit-state \$> lq (c)) is realized in the experiment 
as a physical state | \E f2<J ' (c) ) of N bosons or fermions that 
each carry a (i-dimensional degree of freedom. Since the 
carriers of entanglement are indistinguishable particles, 
each subsystem is of the same dimension d. This degree of 
freedom will typically be an internal property of the par- 
ticle, such as polarization or orbital angular momentum. 
In the following, we will always refer to entanglement be- 
tween these internal degrees of freedom when speaking of 
entanglement between particles. Additionally, the parti- 
cles are distinguished by some external, discriminating 
degree of freedom - typically the spatial mode. Our de- 
scription thus requires N spatial modes (as many as par- 
ticles), such that the minimal size of the single-particle 
Hilbert space for a description in second quantization has 
dimension n = d ■ N, since each particle can live in any 
of the N modes, due to indistinguishability. In order to 
incorporate setups that treat different internal degrees of 
freedom in a distinct way (the internal state of an incom- 
ing particle may change during the scattering process), 
we choose a notation in which the internal (entangled) 
and external (discriminating) degrees of freedom are not 
explicitly distinguished, which will considerably simplify 
the subsequent discussion. For this purpose, the first d 
physical modes (of a total of n modes) are interpreted as 
the d internal states of a particle that is prepared in the 
first spatial mode, the second d modes represent the d 
internal states of a particle in the second spatial mode, 
and so on. In other words, we consider the entanglement 
between the particles found in groups of d modes. 

The creation operator d\ l creates a particle in the fcth 
spatial mode and in the Zth internal state, and we identify 



it - At 



"k,l 



*d(k-l)+V 



(4) 



With this convention, a second-quantized state pos- 
sesses an interpretation in terms of N (first quantized) 
particles that are controlled by N parties - the usual 
quantum information paradigm of an ./V-qudit-state - 
only if there is exactly one particle per group of d modes. 
We call such physical states that can be interpreted as N- 
qudit-state post-selected states, since these are precisely 
the states that are obtained by post-selection. Any N- 
qudit-state in first quantization possesses a representa- 
tion in second quantization, but not vice-versa. 
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The second-quantized version of the initial iV-qudit- 
state |$ l9 (c)) in is thus necessarily a post-selected 
state, 



d 

i* 29 w> = £•••£ ch,.., jN n i vac ) > ( 5 ) 



JV 



j'i=i jV=i 



fe=i 



where |vac) denotes the vacuum. The occupation of the 
n = d ■ N physical modes is characterized by a mode 
occupation list, 



?= (n,r 2 ,...,r„), 



(6) 



which describes a state with rj particles in the jth mode. 
Equivalently, we can specify the modes in which the N 
particles are located, i.e. define a mode assignment list. 
For the state ([5]), the latter reads 

d(r) = (j u d + j 2 ,2d +j 3 ,...,(N- l)d + 3N ). (7) 

We illustrate the mode occupation and mode assignment 
lists of several distinct physical states and the corre- 
sponding iV-qudit-states in Table [I] 



d N 



occupation r assignment d(r) |$ q ) 



3 2 (1,0,0,0,1,0) 

2 3 (QiXdX?) 

2 4 (i^, / m, / m, / ]Co) 

2 3 

3 2 (i7m,oTo7o) 
2 2 (2^0, 0^0) 



(1,5) |1,2> 

(1,3,6) |1,1,2) 

(1,4,6,7) |1,2,2,1) 

(1,3,4) no 

(1,3) no 

(1,1) no 



TABLE I: Physical Fock states and interpretation as 
7V-qudit-states |$ 1<J ). The Fock states are defined by the 
mode occupation list r, Eg. or, equivalently, by the mode 
assignment list d(f), Eq. an. If exactly one particle is present 
in each group of d modes (which are emphasized by the brack- 
ets), there is a corresponding first-quantized iV-qudit-state 
The last three states correspond to unsuccessful scat- 
tering events and cannot be used for quantum-information 
purposes within our scheme. 



As anticipated above, we define an explicit map Q a be- 
tween the first-quantized iV-qudit-state and the second- 
quantized representation, 



d d 



N 

n« 

fc=l 



,t 



| vac) (j 1 ,...,j N \ lg ,(8) 



such that 



|v^(c)) = o |^( c )), 

|a> l9 (c)) = ftt|* a «( c )) = A- 1 |* a «( c )). (9) 



\k=N 



it 

*<i(A!-l)+J t 




d(*-l)+3* 



ao) 



The map f2 a is defined for a given particle-number N, 
and invertible only for the subset of post-selected many- 
particle states. The post-selection procedure is described 
by the projector Pi ja on the subspace of states with ex- 
actly one particle in each group of d modes, 

/ l \ / N 

h,— j'jv = i 

where the inverted order of the indices for the creation 
operators avoids sign changes due to the exchange of 
fermions. For example, the first three states given in 
Table [i] are eigenstates of A, a and possess an interpreta- 
tion as 7V-qudit-state, whereas the last three do not and 
lie in the kernel of P\ a - 



C. Many-particle evolution 

The manipulation of entanglement by the thus defined 
setup relies on the scattering of non-interacting parti- 
cles that carry the possibly entangled internal degree of 
freedom. In other words, the time-evolution of a physi- 
cal initial state |\E' 29 (c)) = ^ Q |$ 1<? (c)) is governed by a 
single-particle scattering matrix W that connects n in- 
put modes to n output modes b\. The state of each 
particle prepared in an input mode thus evolves into a 
superposition of amplitudes localized at different output 
modes: 



if 



E^ 



(11) 



fc=i 



The physical scattering process is governed by the 
Schrodinger equation, it is thus unitary. The scatter- 
ing matrix W , however, is non-unitary, in general, since 
particle loss may occur, i.e. a particle that is prepared in 
an input mode has a finite amplitude to not reach any 
output mode. The matrix W is, in general, embedded in 
a unitary matrix of dimension (2n — 1) x (2n— 1) [60ll61j. 
such that 



U 




(12) 



where X, Y, Z are matrices of dimension n x (n — 1), 
(n — 1) X n and (n — 1) X (n — 1), respectively, i.e. at 
least n — 1 auxiliary modes are needed to account for 
particle loss. The non- vanishing matrix elements of X 
describe the loss of particles, i.e. they quantify the tran- 
sition amplitudes that describe the processes in which 
an injected particle is not collected by any of the first n 
modes. For example, particles may be lost in free-space 
propagation schemes 33, 34, 62 . The embedding matrix 
U can be inverted, and the inverse scattering process on 
the n modes of interest (with the input and output modes 
exchanged) is thus described by W\ and 



*} -> E W IA 
fc=i 



(13) 



fc=i 
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That is to say, even though W does formally not nec- 
essarily possess an inverse, the scattering process with 
exchanged input- and output modes is always described 
by 

The single-particle evolution governed by the ma- 
trix W applies for every single particle. The many- 
particle state \^> 2q (c)) thus evolves according to the time- 
evolution operator M(W), which maps a Fock state of 
particles that were prepared in the input modes, i.e. a 
state given in terms of input-mode creation operators 



a k , to a state in terms of output creation operators &j . It 
thus describes the effect of the single-particle evolution 
W on the many-particle state. Inserting ( 11 ) into (J5J), we 
find the effect of M(W), 



M(W) |* 29 (c)) 



(14) 



N / n 



= J2 c h,~,3N II \ H W d-{k-i)+ Jk .ib\ |vac) 



31=1 3w=l 



k=l \l=l 



For each summand of the outer multi-sum (which runs 
over ji, . . . j N ), we find 



N / n 



(15) 



k=i \i=i 



N 



S S n B/F<» II W d(k-1 



)+3m.o-(m) 



rn—1 



where the sum over f runs over all possibilities, given 
W, to distribute the particles among the final modes, 
d(f) is the mode assignment list, and a runs over all per- 
mutations of this list. The fermionic/bosonic nature of 
the particles is taken into account by sgn B / F (a) , which is 
unity for bosons (B) or the signature of the permutation a 
for fcrmions (F). The time-evolved state M(W) \$> 2q (c)) 
is thus a coherent superposition of all the possibilities to 
distribute the particles among the output modes. Simi- 
larly to single-particle paths that contribute to an event 
that exhibits interference, many-particle paths are super- 
posed coherently [IS]. We illustrate the process with a 
scattering setup for N — 4 particles that each can popu- 
late d = 3 internal states, in Fig. [I] 



1. Properties of the final state 

After scattering through the setup, the final, second- 
quantized state ( 14 ) consists of a coherent superposition 



of all particle arrangements r, with corresponding am- 
plitude given by ( |15[ ). In general, not all components 
can be interpreted as AT-qudit-states, since double oc- 
cupation of some mode or some group of d modes may 
occur. We thus need to project the state onto the sub- 
space of post-selected states via the application of -Pi,6, 
and subsequently apply the inverse mapping Q7 1 , where 



-N f 

' 1 ,2 a l,3 a 2,l 



~^ f — 



x 2,2 "2,3 "3,1 "3,2 "3,3 "4,1 "4,2 



At 



I'll 1**1 | 



w 



^1,1 ^1,2 &1,3 ^2,1 ^2,2 ^2,3 ^3,1 ^3,2 ^3,3 ^4,1 ^4,2 ^4,3 



FIG. 1: (color online) Scattering setup for N = 4 qutrit- 
like (d — 3) particles. The input state corresponds to 
f = (1,0,0,0,1,0,0,0,1,1,0,0), d(r) = (1,5,9,10), and 
|^(c)) = 11,2,3,1). 



fib is defined in full analogy to fl a in &8l, but with inter- 
changed input- and output-creation operators a, and bl . 
We eventually obtain the final iV-qudit-state 



l^Cg)) = Cl7 l P hb M(W)n a \^(c)) 



(16) 



with a norm that equals the success probability of the 
procedure (in the following, unnormalizcd tensors will 
be denoted as such by a tilde). The normalized state 
\$ lq (g)) describes a d-dimensional, iV-partite state, and 
it is fully defined by its coefficient tensor gj 1 ...j N in the 
representation analogous to ([I]), which allows us to char- 
acterize its entanglement properties with the tools of 
quantum information theory [171 158] . 



D. Entanglement creation and detection 

The above-defined mapping between first- and second- 
quantized representation and the general procedure for 
the manipulation of entanglement via many-particle scat- 
tering are illustrated in Fig. [2] As we show now, the 
process can be used as a protocol for entanglement gen- 
eration, detection and projection: 

Entanglement generation. Given an initial separable 
state |$ l9 (c)), a target entangled state 



\$ lq (g)) = ft" 1 P hb M(w)n a \^(c)) 



(17) 



is generated by scattering and post-selection. 
The success probability of the entanglement gen- 
eration process is given by (& lq {g)\& lq (g)) = 

P lfi M(W)(l a |$ l9 (c)) . For example, Refs. 0[371 
l3"8] rely on this procedure. 

Entanglement detection. Given an entangled ini- 
tial state |$ 1? (c)), the scattering process yields 
a superposition that contains one or sev- 
eral separable (known) signal states \^ lq (s^)), 
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FIG. 2: (color online) Entanglement manipulation 
via many-particle scattering. Given the iV-qudit-state 
|$ llJ (c)}, Eq. nTb, the physical second-quantized state is given 
by Eq. jgj, |*^(c)} = D. a |$ 1<J (c)), with the map operator 
fi a given by Eq. dsb - The scattering of t he p articles through 
the setup is described by M(W), Eq. (14 1, the projection 



on the subspace of post-selected states is mediated by Pi,b, 
Eq. (10 1. The mapping to a first-quantized state occurs via 
Cl' 1 , Eq. The sketch on the right-hand side shows the 
physical situation for d = 4 and N = 2. The initial state (up- 
per sketch) contains only one component, it is injected into 
the scattering setup. The final state contains many coherently 
superposed components, some of which are not interpretable 
as iV-qudit-states (marked with red crosses). The projector 
Pi,b describes the post-selection mechanism, i.e. it eliminates 
the unwanted components and leaves only the post-selected 
components intact (marked by green crotchets). 



i.e. \<S> lq {~g)) = fii'A.lM^fi.l^'fc)) with 
{&<i(sW)\® lq (g)) ^ 0. All states that are orthog- 
onal to the initial state do not lead to any signal 
state, such that the detection of signal states is a 
signature for the initial entangled state. 

Entanglement projection. If we start with a superpo- 
sition of a desired target entangled state \§ lq (g)) 
- as for the entanglement detection scheme - with 
other undesired components, the detection of a sep- 
arable signal state \^ lq (s^)} projects the initial 
superposition onto the target state. This proce- 
dure is used, for example, in entanglement swap- 



ping protocols 23!, or to project atoms into entan- 
gled states by measuring the photons emitted by 
the atoms [53] , 

The three tasks pose very distinct experimental re- 
quirements in general, and they are therefore typically 
not realized in the same experimental setup. Still, they 
are fully equivalent in theory: For any number of particles 
iV and any dimensionality of the entanglement carrier d, 
every created state can also be detected, and every state 
that can be detected can also be created, in principle. 

More formally, a target entangled state \$> lq (g)) that 
is created via the scattering of an initial, separable state 
|<fr l9 (c)) on a setup W can be detected with the very same 
separable state as the signal state, via the scattering of 
the target entangled state on W\ i.e. using the output 
modes of W as input. 

In order to show this, we need to ensure that \^> 2q {c)) 
qualifies as a signal state. We thus require (i) a finite 
detection probability, 



(^ 2q (c)\M(W^)\^ 2q (g)) t^O. 



(18) 



Furthermore, (ii) for all states characterized by g that 
are orthogonal to the target state g, the signal state is 
never produced, i.e. 



V 3 \ with (<5> lq (g ± )\<5> lq (g))=Q: 
(y 2q (c)\M(Wi)\V 2q (g ± )) = 0. 



(19) 



Eqs. ( 18 ) and ( |T9[ ) are actually satisfied as a direct con 
sequence of M(W r ) = M*(W), i.e. by the fact that the 



inverse scattering process is described by W* (see ( 12 ) 



and (|13j)). 

In fu 1 analogy, any state that can be detected can 
also be created. In order to facilitate our notation in the 
following, we will always refer to the entanglement gener- 
ation of a desired target entangled state - our results can 
then be applied analogously for entanglement detection 
and projection. 



E. Treatment of internal degrees of freedom 

The principle mode of operation being established, it 
remains to specify the properties of the scattering matrix 
W for specific physical systems. 

The formalism can accommodate setups in which par- 
ticles are simply redistributed among output modes, 
without any manipulation of their internal degrees of 
freedom. We will refer to such setups as non-polarizing 
setups (which does not necessarily imply the restriction 
to polarization as the internal degree of freedom). Se- 
tups that are sensitive to the internal degree of freedom, 
i.e. polarizing setups, are also thinkable (such as polar- 
izing beam-splitters, which redirect photons conditioned 
on their polarization, but which do not change the po- 
larization of the incoming particles). We call a setup 
that also manipulates the polarization, e.g. by using wave 
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plates, polarization-manipulating setup. For these three 
different types of setups, the scattering matrix W will as- 
sume a distinct structure, and distinct target states are 
accessible. 



1. Non-polarizing setups 

When photons propagate in free space or scatter off a 
non-polarizing multiport beam-splitter that treats all in- 
ternal degrees of freedom equally, we have a propagation 
of the form 



N 
k=l 



(20) 



for particles that are prepared in a spatial mode j and in 
the internal degree of freedom I, where we distinguished 
internal and external degrees of freedom for the moment 
(See Eq. Q). Each de gree of freedom is treated inde- 
pendently and equally, which formally translates into a 
scattering matrix W that contains one elementary scat- 
tering matrix V of dimension N x N, 



W 



V 



(21) 



For example, the matrix that describes a non-polarizing 
beam-splitter, 



W 










1 


o ^ 




1 





1 





1 


1 




1 





-1 











1 











® 1(>22) 



possesses precisely the structure (21 1. Any scheme that 
does not manipulate the internal degree of freedom of 
the scattered particles, such as the ones discussed in [33T 
1551 1571 155] , can be accomodated by a matrix of the form 



The entanglement content of the final TV-qudit-state 
depends on the prepared initial internal states of the par- 
ticles: The multiset of internal states of the particles is 
not altered, i.e. an initial state = |ei, e%, . . . ejv) will 
lead to a final state that can be expressed by superposi- 
tions of terms which consist of permutations of the |e^) on 
the output modes. If all \ek) are equal, no entanglement 
can be created. 



2. Polarizing setups 

The situation is similar for polarizing setups, where 
we allow different dynamics for each internal state, but 
no coupling between them, e.g., a horizontally polarized 
photon may be directed to another output mode than a 
vertically polarized photon, but the polarization of in- 
coming photons is never changed. This model is real- 
istic for degrees of freedom that can be discriminated, 



but which cannot be manipulated at will, i.e. the device 
equivalent of a wave plate is not available. The scattering 
matrix is still a direct sum of scattering matrices, which, 
however, do not need to be equal, 

W = F (1) ®diag[(l,0...,0)] 

+V {2) ® diag[(0, 1, . . . , 0)] 



+V M ®diag[(0,...0,l)], 



(23) 



where the are N x N matrices and diag[v] is the 
matrix that contains v on the diagonal, whereas all other 
elements vanish. 

For example, a polarizing beam-splitter does not mod- 
ify polarization, but its action on the modes depends on 
the polarization state. The matrix 



Wpbs 




(24) 



®diag[(l,0)] 



diag[(0, 1) 



describes a polarizing beam-splitter which reflects and 
transmits vertically polarized photons with the same 
probability and always transmits horizontally polarized 
photons. 



3. Polarization-manipulating setups 

Finally, a setup which does not only treat different 
internal degrees of freedom in a distinct way but also 
allows their manipulation is described by an unrestricted 
scattering matrix W. This is experimentally more de- 
manding in general, but routinely done for photon polar- 
ization, with the help of phase shifters, polarization ro- 
tators (e.g. implemented by the combination of quarter- 
and half-wave-plates) and polarizing beam-splitters. 



F. Effective degrees of freedom and independence 
of initial state 

1. Multiple mode occupation 

At first sight, the formalism seems to restrict severely 
the types of setups with which entanglement is produced: 
Often, the initial state in an experiment does not possess 
an interpretation as ./V-qudit-statc in the quantum infor- 
mation abstraction, such that no |$ 1<J ) can be given. For 
example, a many-particle state with all particles in the 



same spatial mode is scattered on several beam-splitters 
in |14j . which permits the creation of a permutation- 
symmetric entangled state. 

Such initial multiple population is, however, unprob- 
lcmatic: If two particles are prepared in the same mode, 
we can equally well describe them by a setup in which 
they formally occupy two distinct modes, but for which 
the two rows of the matrix representation that describes 
the evolution of the particles are identical. We can thus 
absorb such initial preparation of particles into the ma- 
trix W. 



2. Local unitary operations 

All local operations that act independently on each 
particle in the input or in the output modes, i.e. local 
unitary operations [53J |M] , can also be absorbed by the 
scattering matrix: 



W = 



(4° 



7 ° (4 



(/) 



L 



(25) 



where and are d x e?-matrices that describe local 
unitary transformations on the fcth spatial input- and 
on the Zth spatial output mode, respectively. We can 
therefore restrict our analysis to a unique initial state 
= |1, 1, 1, . . . , 1), in which all particles are prepared 
in the same internal state and in distinct spatial modes. 
All relevant parameters can be varied by the entries of 
the scattering matrix W . 



III. STATE REPRESENTATION AND 
COMBINATORIAL BOUNDS 

Having established the physical process for entangle- 
ment generation as well as the possible structures of the 
scattering matrix W , we now proceed to the explicit rep- 
resentation of the states that can be generated, given 
an initially separable state. The form for polarization- 
manipulating setups will be the most general one, while 
we find a simpler representation for non-polarizing se- 
tups. In all cases, elementary algebraic manipulations al- 
low us to obtain a form from which the maximal Schmidt 
rank, i.e. the maximal number of separable components 
of a state, can be read off. 



Non-polarizing setups and permutation 
representation 



where |efc) denotes the internal state of the fcth particle. 

Since, by assumption, the internal degrees of freedom 
cannot be modified by the setup, the multiset of inter- 
nal states that the particles are initially prepared in, 
{ | ei) , . . . |e./v)}, remains invariant under scattering. Con- 
sequently, the first-quantized final state will always pos- 
sess a permutation representation. 



\& q {9)) 



£ a. 



(27) 



where g a is a coefficient that depends on the permuta- 
tion it, and all possible distributions of the single-particle 
states \ej) on the spatial output modes are taken into ac- 
count. By a parameter-counting argument, we see that 
not every state can be represented by the form (27 1: The 



number of parameters (the number of g a ) is N\, whereas 
N N parameters are required for the description of a gen- 
eral state of N particles each associated with an N- 
dimensional degree of freedom. For constant coefficients 
9<t — 9> the state (27 1 is fully permutation-symmetric. In 
particular, any permutation-symmetric state with d = 2 
can be written in the Majorana-representation |651 166) . 
i.e. one can then always find a multiset {|ej)} to rep- 
resent it in the form (27). Physically speaking, every 



permutation-symmetric state of qubits can be created 
with beam-splitters and phase-shifters (36| . 

In the present scattering scenario, the coefficients g a 
cannot be adjusted independently, but they are given as 
a function of the scattering matrix elements Vk t i, and 
they depend on the type of particles (bosons or fcrmions) . 
Each component le^m, . . . , € a (N)) m the state (27l ob- 
tains the amplitude of the corresponding many-particle 
path, for which the cr(j)th particle falls into the jth spa- 
tial output mode. Consequently, applying to the 



state in (15 1 and identifying the coefficients g a in the 



equivalent representation (27) yields 



N 



g a = sgn B/F 



(28) 



where sgn B y F (cr) denotes the signature of the permuta- 
tion a for fermions, and is always unity in the case of 
bosons. That is to say, the final state for non-polarizing 
setups is a superposition of all possibilities to distribute 
the particles prepared at the N spatial input modes over 
the N spatial output modes. The amplitude for each 
of these permutations is given by the product of the 
scattering matrix elements as defined in ( 28 1 . 



For a non-polarizing setup, we are given a scattering 
matrix of the form (21), with an N x A- matrix V as 



non-trivial ingredient. Any separable initial state can be 
written as 



(26) 



B. Polarizing and polarization-manipulating setups 

For polarizing setups, which do not manipulate the 
internal degrees of freedom of the particles, one cannot 
use the permutation representation to express the final 
state, since particles in distinct internal states can behave 
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(a) 



(b) 



w 

1 1 1 1 1 1 1 
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For the interpretation of (30), we can apply an intu- 
itive picture of many-particle paths [48], which is also 
illustrated in Fig. [3j Many distinct many-particle paths 
contribute to post-selected events with one particle per 
group of modes, and the entanglement in the final state 
originates from this superposition of distinct paths [57] 
(see Fig. [3] (a)). Each separable basis state \ji,j2, ■ ■ -Jn) 
in the final state representation is fed by up to AM dis- 
tinct contributions (see (b) and (c)): The particle ini- 
tially in mode d ■ (a[k) — 1) + 1 is redirected to mode 
d-(k— l)+jk, and the product of the corresponding ampli- 
tudes gives the many-particle amplitude that contributes 
to this specific scattering process. All these amplitudes 
are superimposed and interfere (see (b)), and their sum 
equals the coefficient gj lt ...j N - 



FIG. 3: (color online) Superposition of many-particle 
paths, (a) A state of N = 2 particles that carry a d — 4- 
dimensional degree of freedom is characterized by the coef- 
ficient tensor g, which indicates the amplitude of each basis 
state, (b) Each basis state (here: |1, 1}) is fed by the many- 
particle paths that contribute to the respective event, (c) In 
general, there are up to AM of these paths, we show here two 
exemplary paths for N — 4 particles. 



differently and the single-particle states {|efc)} are not 
preserved. We therefore treat this case together with 
polarization-manipulating setups. In general, the target 
state's first-quantized notation is 

d d 

l* l9 (.9)> = E ~9h,-,m\h,h,---3N).(29) 

jl = l JN=1 

For an initial state with all particles in the internal state 
\jk = 1) (remember that arbitrary local operations on the 
initial qudits can be absorbed in the scattering matrix W, 
as described in Section [II F| we evaluate the coefficients 
9ji....,j N via the application of Pi j, onto (15): 



9h,...,j N = (30) 

TV 

SglES/F^) II W / d-( ( r[fc]-l) + l,d(fc-l)+ J)! - 



E 



k=l 



C. Combinatorial bounds 

With the explicit state representations at hand (see 
(27) and (30)), we can now map the task "given a cer- 
tain infrastructure, can we create the state |$ lg (/i)) start- 
ing from a separable state?" to a concrete mathematical 
problem: It corresponds to the set of d N equations 



9ji. 



(31) 



for the matrix elements Wkj, where the gj t ,...j N are given 



by ( 30 1 and r\ accounts for the sub-normalization of the 
g. In general, this task is manifestly rather difficult, since 
the coefficients 9j 1 ...j N are highly non-linear functions of 
the unknown parameters Wk,i- 



Still, the explicit representation (30) permits the iden 



tification of limitations that apply to entanglement gener- 
ation via many-particle interference, in terms of an upper 
bound for the generalized Schmidt rank [49j . 



The bound can be derived from the representation of the final state \<fe lq (g)), which reads, according to ( 29|30 l, 

d d N 

\® lq (g)) = E " E E S & n B/F ( CT ) II W d(<r[k]-l)+l,d{k-l)+j h lil) ■ • • Jn) , (32) 



ii=i j«=icreS]v 
and which can be re-written as 



k=l 



l* l9 (5)) = E S S n B/F( cr ) 
<t€Sn 



1 E Wd(alm]-l) + l,d{rn-l)+j m \jm), 



(33) 



where \j m ) rn denotes the state of the 771th qudit prepared in the j m th internal state. In other words, every state that 
is created by many-particle interference can be written as the sum of AM (from the sum over all permutations a) 
separable terms of the form 
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This result immediately yields an upper bound for the 
Schmidt rank R(\$ lq (g))) of the state \<$> lq (g)), 



R(\<S> lq (~g)))<N\ 



(34) 



where the Schmidt rank is the minimal necessary num- 
ber of summands in the representation of \Q lq ) as a 
sum of separable states [IS]. This fundamental bound 
for multipartite entangled states that are generated via 
many-particle interference can be understood via an in- 
tuitive combinatorial argument: Each possibility to dis- 
tribute the particles among the output modes can only 
contribute one term in the representation of the emerg- 
ing state. Thereby, we have established a link between 
the number of many-particle paths and an entanglement 
measure, namely the generalized Schmidt number, which 
is defined as \og 2 R(\$ lq (g))) gS_]. 



D. Initial states with multiple occupation 



Through our reference to ( 30 1 , we assumed in Sec- 
tion |III C| above that the initial state is the iV-qudit- 
state = |1,...,1), i.e. a state with exactly one 

particle in each group of d modes. If an initial physi- 
cal state with larger occupancy is used (which does not 
admit an interpretation as iV-qudit-state and which is 
only admissible for bosons), a tighter bound emerges. 
In this case, we can formally use the same initial state 
but we account for the multiple occupation by 
equating the corresponding rows of the scattering matrix 
(see Section II F ) . When rj is the occupation number of 
each mode (zJ7 = i r j = N), the scattering matrix W ef- 
fectively reflects this multiple occupation by possessing 



rj equal rows. The outer sum in ( 33 1 reduces to only 
■An/(n7=i r J'0 terms, since for each a £ Sjy, there are 
n^ =1 rj\ permutations that leave the summand 



/ •• d(tr[m]—l)+l,d{m—l)+j m \3m) 
3m=l 



(35) 



invariant. This can, again, be interpreted in a picture of 
many-particle paths: When several particles are initially 
prepared in the very same state, a permutation of these 
particles does not lead to a distinct many-particle am- 
plitude. Therefore, the number of distinct many-particle 
amplitudes is reduced by a factor J^J - =1 rjl. 

As an extreme example, it is not possible to extract 
entanglement from an ideal Bose-Einstein-Condensate 
(BEC), which is described by 



I* 



2q v _ 
BEC/ — 



(< 



vac) 



(36) 



i.e. a state of N particles with all particles in the same 
mode, since N\/\]\ n j=l rj\) = 1 for f = (N, 0, . . . , 0). 
Post-selected events with one particle in each output 



mode can only be reached via one many-particle path, 
such that the Schmidt rank of the resulting state is al- 
ways unity. This generalizes the result of [67], where it 
was found that no entanglement can be extracted from 
a BEC when no manipulations of the internal state dy- 
namics are allowed. One consequence of our results is 
that even such manipulation on the internal degrees of 
freedom is of no help. 



E. Entangled initial states 

In turn, the use of an entangled initial state allows 
one to increase the upper bound on the Schmidt rank of 
the final state: The combinatorial bound applies to every 
single separable term in the representation of the initial 
state in its generalized Schmidt form. Given an initial 
Schmidt rank R m , the combinatorial bound is extended 
up to R m N\. In other words, the concatenation of many- 
particle interference processes and post-selection can be 
used to increase the Schmidt rank step by step-wise, by 
a factor of up to Nl in each step. 



F. Interpretation 

The combinatorial bound should not be regarded as a 
purely quantum phenomenon: It also applies to a classi- 
cal probabilistic process in which N distinguishable ob- 
jects are randomly distributed to N observers. The num- 
ber of distinct outcomes that the N observers can expe- 
rience is, again, bound by Nl, just like the generalized 
Schmidt rank is according to (34 1. If some of the objects 
are identical, permutations that exchange these objects 



do not lead to inequivalent events, just like the multi- 
ple initial occupation of any mode reduces the rank, as 
discussed in Section |IIID| The quantum properties of 
bosons and fermions, i.e. their (anti)-commutation rela- 
tions which provide an alternating sign for fermions in 
the sums ( 15|30 ), do not play any role, and the bound 
applies equally well to both species. It is, furthermore, 
independent of the presence of polarizing components. 
In contrast, we will see that bosons and fermions indeed 
do hold a different potential for entanglement generation, 
as will be shown by means of the dimensionality of the 
manifold of accessible states in the next Section. 



IV. DIMENSION OF THE MANIFOLD OF 
ACCESSIBLE STATES 

For systems of large dimensions d > N, the combina- 
torial bound (34 1 dominates the set of accessible states. 



By a parameter-counting argument, however, we see that 
not all states can be created when the number of parti- 
cles N is large, for any dimension d > 2: On the one 
hand, the number of complex parameters for a general 
iV-partite quantum state with a e?-dimensional internal 
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degree of freedom is d N (including normalization and 
global phase). On the other hand, the number of com- 
plex paramete rs of the scattering matrix W that controls 
the final state (32) is only d-N 2 : We are given the initial 
state \<fr lq ) = |1, 1, . . . , 1), and only the matrix elements 
W(k-i)d+i,j w ith 1 < k < A and 1 < j < dN appear 



in the sum (32 1. The exponential growth of the param- 
eters of entangled states with the number of involved 
particles reflects the complexity of many-particle quan- 
tum mechanics, which cannot be reproduced by the lim- 
ited amount of free parameters in the interference-based 
setup under investigation. 



A. Fermionic and bosonic bounds 

The rough upper bound provided by the above 
parameter-counting argument can be ameliorated by 
studying the structure of the state coefficients in (32) 
in more detail. This allows us to obtain a tight upper 



bound for the dimensionality of the manifold of accessi- 
ble states. We denote by S F ( B ) the manifold of states of 



the form ( 32 ) for fermions (bosons) , which is spanned by 
all matrices W. In contrast to the combinatorial bound, 
we find a striking difference between bosons and fermions. 
The derivation of the bounds is presented in Appendix 
IbI For fermions, we find 



dim(S F ) < (d - 1) A 2 - A + 2, 
while for bosons, we have 

dim(S B ) < dN 2 - 2 A + 2. 



(37) 



(38) 



Roughly speaking, the computation of scattering ampli- 
tudes for fermions involves the evaluation of determi- 
nants, which is facilitated by the symmetry properties of 
this function, in particular, by the property det(AB) = 
det(A)det(_B). Scattering amplitudes for bosons instead 
involve the permanent |68H70| , a complex- valued function 
of matrices similar to the determinant, but that does not 
feature equivalent symmetry properties. The number of 
effective parameters for bosons is therefore larger than 
for fermions. 



The bosonic bound (38) is only applicable when it is 
smaller or equal to the dimension d N of the total state 
space. For N = 2, and for d = 2, A = 3,4, the bound 
exceeds the total dimensionality of the many-particle 
Hilbert-space, N d ; it is thus not relevant. For large sys- 
tems, the bound strongly limits the manifold of accessible 
states, as we will also see in the next Section in Table |H| 



B. Tightness of bounds 



The established inequalities (37) and (38) provide up- 
per bounds to the dimensionality of the manifold of 
entangled states that can be created by fermions and 
bosons, respectively. For small dimensions d and small 



particle numbers A, we found numerically that the equal- 
ity holds in (37) and (38). One thus indeed possesses a 



number of independent parameters that corresponds to 
the right-hand side of (37), (38). For this purpose, we 



evaluated the rank of the Jacobian 



(39) 



for random scattering matrices W. The results are shown 
in Table [XTJ they always match the encountered upper 
bound. 





TV Bosons 






TV Fermions 




2 3 4 5 


6 7 


8 


2 3 4 5 6 7 


d = 2 


4 8 16 32 


62 86 


114 


4 8 14 22 32 44 


3 


(8) 23 42 67 


98 




8 17 30 47 68 


4 


(12) 32 58 92 


134 




12 26 46 72 


5 


(16) 41 74 117 






16 35 62 97 



TABLE II: (color online) Numerically obtained lower 
bounds for the dimensionality dim(EB/F) of the acces- 
sible set of states, for unconstrained setups, according 

to Eq. ( 39 ). Blue underlined numbers indicate combinations 
of d and N which are not subject to any constraint, i. e. all en- 
tangled states can be created via many-particle interference. 
Red numbers (in parenthesis) refer t o co mbinations of (d, N) 
for which the combinatorial bound (341 is tighter than the 
bosonic bound ( 38 1 . All numerically inferred lower bounds 
match the analytically obtained upper bounds (1371) and 1381) 



or the maximal dimensionality that can be attained according 
to the combinatorial bound (341. 



C. Interpretation and consequences 

Our results on the effective number of parameters for 
bosons and fermions demonstrate that the combinatorial 



bound (34) is not the only restriction to the set of acces- 
sible states. The complexity of general multi-partite en- 
tangled states cannot be reproduced in interference-based 
setups, simply due to the limited number of parameters. 
When fermions are used, the symmetry properties of the 
determinant reduce the effective number of parameters 
in comparison to bosons. Indeed, the computation of 
fermionic many-particle scattering amplitudes is compu- 
tationally much less demanding than the corresponding 
task for bosons and this computational complexity 
is reflected in the larger manifold for bosons. 

This larger manifold of states that are accessible for 
bosons comes at an expense: For fermions, the entangle- 
ment generation process is assisted by the Pauli princi- 
ple, which forbids multiply occupied modes that corre- 
spond to unsuccessful events. The average tendency of 
bosons to bunch, on the other hand, enhances unsuccess- 
ful events. For larger dimensions d, the total number of 
modes n increases, and the Pauli principle becomes less 
dominating, as well as the bunching tendency of bosons is 
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less prominent: Due to the lower particle density, multi- 
ply occupied states become less probable in general, since 
many more possible states with at most one particle in 
each of the n modes emerge. Simultaneously, the quo- 
tient of the bounds to the manifold dimension of fermions 
to the manifold dimension of bosons, dim(Sp)/dim(SB) 
(see (37) and (38)), converges to unity for d —> oo, i.e. the 



advantage of bosons with respect to fermions regarding 
the dimensionality is jeopardized. In other words, there 
is a trade-off between success probability and the dimen- 
sionality of the accessible manifold. 



Given the dimensionality estimate ( 38 1 for qubits that 



are carried by bosons, we can expect that all four-partite 
and five-partite qubit states can be created with linear 
optics, with experiments similar to [13j HH [18j [31] [71]. 
For six photons, we find that only 62 of a total of 64 
dimensions are spanned. Given the manifold possibili- 
ties to create entangled six-photon states by parametric 
down-conversion processes, it is, however, possible that 
one may find ways to create the unaccessible states start- 
ing from entangled initial states. 

Significant progress has been achieved in the design of 
alternative, efficient single-photon sources [72TfT4"] . How- 
ever, our results indicate that extending experiments to 
larger particle numbers is not only limited by technolog- 
ical issues, but fundamental limitations to the accessi- 
ble states will severely come into play for large particle- 
numbers. For example, for N — 8,d — 2, Sb has only 
114 dimensions, significantly less than the dimension of 
the HUbert space, 256. In analogous experiments with 
fermions, e.g. with cold atoms in optical lattices, one will 
experience a bound for the accessible states already for 
four particles. 

For bosons, the inclusion of loss in the setup, i.e. the 
use of non-unitary scattering matrices W , permits the 
generation of states that are inaccessible for unitary ma- 
trices: Consider a lossless setup described by the unitary 
dN x dN matrix W, i.e. the matrices X, Y, Z in the rep- 



resentation ( 12 ) are empty. The number of parameters 



that affect the relevant dN x iV-dimensional sub-matrix 
W' of W can be inferred by counting the Euler-angles 
of W that do not affect W'. The total number of com- 
plex Euler-angles that define W amounts to Nd(Nd — 1), 
while N(d — l)(N(d — 1) — 1) of these are irrelevant for 
W' . The difference yields the number of parameters K, 



A- = u-iV-» 



(40) 



This value is exceeded by the dimensionality of bosons, 
dim(Se) (38). The inclusion of lossy setups (or, equiva- 



lently, of the multiple population of modes in the initial 
state, see Section IIID[ ) thus offers advantages over loss- 
less, unitary setups. For fermions, however, the number 
of parameters of a unitary matrix U is larger than our 
dimensional bound, dim(U) > dim(SB), which suggests 
that no advantage is accomplished for fermions by in- 
cluding loss in the setup. 



V. CONCLUSIONS AND OUTLOOK 

The emergence of quantum correlations from the prop- 
agation and detection of bosons and fermions can be de- 
scribed by the presented unified formalism, which also 
provides an intuitive interpretation: Particles are re- 
distributed from input modes to output modes, and the 
probability to fall into each output mode depends on the 
initial state of the particle. This process has a purely clas- 
sical analogy: The random distribution of distinguishable 
objects to parties/observers also leads to classical correla- 
tions between the objects that are subsequently observed 
by the parties. Here, we have considered the coherent 
version of this process, and the resulting correlations are 
of quantum nature. 

The physical process of redistributing particles among 
observing parties is the very basis for a panoply of exper- 
iments and theoretical schemes. We aimed at providing a 
common framework for the creation and detection of en- 
tangled states via the deletion of which-way information, 
and obtained an explicit representation ( 30 ) of all states 



that can be created or detected via linear optics setups. 
Depending on the experimental capabilities, i.e. whether 
the device that is equivalent to a polarizing beam-splitter 
can be implemented or not, the form of the scattering 
matrix is very different, and different types of states can 
become accessible. The combinatorial bound, however, 
cannot be circumvented. The computational complexity 
of the permanent and determinant, which appear in the 
scattering amplitudes for bosons and fermions, respec- 
tively, is reflected by the effective number of parameters 
that are available to generate entangled states. 

In general, any experiment in which particle-carried 
entanglement is created without direct interaction be- 
tween the constituents can be described by a scattering 
matrix W, and the resulting state has a representation of 
the form (30). Given the experimental limitations, one 



can thus optimize the creation of entanglement by finding 
the optimal matrix W. 

Some proposals that encode qudits in a different way 
than we have assumed here can circumvent the combi- 
natorial bound, and can reach states that are not repre- 
sented by ( 30 1 . For example, mode-entanglement [75H77] , 
i.e. entanglement in the particle-number degree of free- 
dom, is not subject to the bound, as well as schemes in 
which qudits are carried by several particles instead of a 
single one (2BESHED]. 

The prospect for many-particle entanglement creation 
with photons is restricted by the dramatic decrease of 
the success rate for increasing particle number N, since 
unsuccessful events with multiply occupied spatial out- 
put modes then become more and more probable, al- 
though there are approaches to tackle this inefficiency 
by using the time-bin as an additional discriminating de- 
gree of freedom [81 j . The very understanding of multi- 
partite entanglement being a widely open issue, the state 
representation ( 30 ) offers a way to understand at least 



the structure of the states that can be created via the 
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exploitation of the indistinguishability of the particles. 
In particular, completely anti-symmetric and completely 
symmetric states are contained in the general states for 
fermions and bosons, respectively. Since these special 
classes are rather well understood, the hope is fed that 
one may characterize the more general state ( 30 ) via 



where (U 



FoiA 



entanglement measures or via its non-local properties, 
and thereby also characterize the potential of non-linear 
single-photon schemes [32] ■ 

In particular, the question arises whether the mani- 
fold of accessible states for fermions is a submanifold of 
the accessible states for bosons. This appe ars sug ges- 
tive from the clear hierarchy of the bounds ( 37|38 1 for 
the dimensionalities that we have encountered, but, in 
turn, also questionable, since fully symmetric and fully 
anti-symmetric states which are accessible for bosons and 
fermions, respectively, when the scattering matrix V is 
the constant matrix (see Eq. ( |A3[ )) possess very distinct 
entanglement properties [83] . Also the influence of re- 
strictions on the setup, i.e. the use of non-polarizing se- 
tups instead of the unconstrained form, remains to be 
studied. The answers to these questions may contribute 
further to an understanding of the structure of multipar- 
tite entangled states. 
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Appendix A: Accommodation of existing schemes 

Many existing schemes that seem rather unrelated to 
each other at first sight can be put onto a common foot- 
ing. In the following, we give an exemplary overview over 
various methods that can be described by the formalism 
described in this article, which stresses the shared un- 
derlying physical roots of the different theoretical and 
experimental methods. 



1. Theoretical proposals 

a. Scattering on the Fourier multipart 

Several multipartite entangled states can be created 
with Fourier multiport devices via the scattering of 
bosons that are prepared in certain initial states |38j . 
The setup is non-polarizing, with 



o^Sd i) } s the Fourier matrix. 



No devices that are equivalent to polarizers need to be 
used, such that the infrastructure is independent of the 
specific degree of freedom that is entangled. The final 
state reads (as given in (27)) 



N 



Fou 



5fe=i \ e <j(k)) k , 



(A2) 



<jGS n \j'=1 



where the |efe) are the initial internal states of the bosons 
at the input modes. In particular, a W-type state [84] 
can be created by preparing one photon in horizontal, 
| if), and the remaining A — 1 photons in vertical polar- 
ization, \V). Due to the cyclic invariance of the Fourier 
matrix, it is ensured that all amplitudes of the distinct 
components with the jth particle in horizontal polariza- 
tion (\H, V, V, V, . . .), \V, H, V, V, . . .), etc.) enter with the 
same amplitude. For four particles, a GHZ-state of four 
photons can equally well be produced with the initial 
state \V,V,H,H). 



b. Free space propagation 

A free-space photon-propagation scheme [53"H3"T] allows 
to generate a variety of entangled states. Both, entan- 
glement projection (with consequent entanglement swap- 
ping) [35] . and entanglement creation [37] are possible. 
In all cases, photons propagate in free space, and they 
are collected by an array of N detectors, each of which 
projects the incoming photon onto a certain polariza- 
tion state. Each photon has the same probability to be 
absorbed by each detector. In particular, Dicke states, 
i.e. states \S,m) which are simultaneous eigenstates of 
the total squared spin operator S* 2 and its z-component 
S z , with eigenvalues S and m, respectively, are accessible 
|33j . More generally, all permutation-symmetric states 
can be created |37j , and - with few modifications of the 
scattering matrix - all 2 N distinct total angular momen- 
tum eigenstates [MJ [83] . 

The scattering matrix W is non-polarizing, and its sub- 
matrix V (see ( plj )) is the constant matrix with Vj t k — 
y/p, where p is the probability that a photon is registered 
by the detector. For the state representation, we can use 



(27) to find 



io - e n^w, 

<res N \j=l 



/ I "3fc=l |e<r(fc)) fc 



-<r(fc)/ 



(A3) 



W = V <■ 



and V = U 



Fou 



(Al) 



In other words, a fully permutation-symmetric state 
emerges 03 H7J EH HEl M\ ■ 
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2. Experimental implementations 

Many qualitatively distinct classes of entangled states 
HI EJ [88] have been realized in experiments that em- 
ploy photons as carriers of entanglement. In order to ex- 
emplify the versatility of many-photon entanglement gen- 
eration, we discuss in more detail the projection onto a 
tripartite entangled GHZ-state [89] and the versatile tun- 
ability of distinct SLOCC classes for four-photon states 



via a GHZ-state analyzer [5U] (dashed area in Fig. EJ). 
The polarizing beam-splitters (PBS) reflect vertically 
polarized photons and transmit horizontally polarized 
ones, while the polarizers after the PBS only transmit 
the |+) = (\H) + |y))/ v / 2-component of the light field. 
The resulting scattering matrix Wghz that describes the 
propagation of the photons in the input modes 

|if>i,l^>i, \H) 3 ,\V) 3 ,\H) 5 ,\V) 5 , (A7) 
to the detectors D\, D3, D5 reads 



a. GHZ entanglement swapping 

Entanglement swapping among three parties was ex- 
perimentally demonstrated |89j with the setup repro- 
duced in Fig. [4] A laser pulse that induces parametric 
down-conversion on three aligned crystals (BBO in the 
Figure) results in the creation of 



|4> J 



(A4) 



'3,4»I T /5,6) 

where |$ + ) fe j denotes the |$ + ) Bell-state in the polar- 
ization degree of freedom, shared between the photons in 
the spatial modes k and I: 



1 



H,H) k , + \V,V) 



Photons 1, 3 and 5 are then projected onto a GHZ-state, 

GHZ state analyzer 




UV-pulse B p°^l-/ 



PBS1 1^ 
BBO 3 





PBS 



FIG. 4: (color online) GHZ entanglement projec- 
tion. The photon pairs (1,2), (3,4) and (5,6) are initially 
polarization-entangled, as described by Eq. ( |A4[ ). By the pro- 
jection of photons 1, 3 and 5 onto the GHZ-state (at detectors 
Dl, D3, D5, by the implementation of the scattering matrix 
(A8l), photons 2, 4, and 6 are also left in the GHZ-state. 
Courtesy of C.-Y. Lu, T. Yang and J.-W. Pan [89]. 



\GBZ)l 



3,5 



V2 



\H,H,H) + \V,V,V)), (A6) 



IFghz = 



V2 



/0 1 \ 

1 

1 

1 

1 

v o o i y 



(A8) 



i.e. it corresponds to a combined action of polarizing 
beam-splitters and polarizers, as described in Section 
HE 2 The state at modes 2,4 and 6 is then measured in 
different bases to verify quantum correlations. The path 
lengths Aril and Ad2 (see Figure) are adjusted by moni- 
toring four-photon a nd si x-photon coincident events. As 
discussed in Section 



II D 



(A5) used for the creation of a GHZ-state. 



the matrix Wqjjz can also be 



b. Family of four-photon states 

The setup for the projection onto the GHZ state dis- 
cussed above projects onto one very specific state. One 
can, however, also scan through entire families of entan- 
gled states by tuning the parameters of the scattering 
matrix W . This was realized in a four-photon experiment 
[3T1 I7T] , in which one starts with a four-photon state cre- 
ated via parametric down-conversion, which reads 



1 

2^3 



it 



"l,H a 2,V 



l l,V a 2,H 



-2a x H fi-L v a 2 H 



3 1 



vac) 



(A9) 



i.e. there are always two photons of horizontal and of 
vertical polarization, respectively, and there is one com- 
ponent with the photons distributed among all four phys- 
ical modes and two components with two photons in the 
same mode. 

A variable polarization-rotation with angle 7 is applied 
to the spatial mode d\ (see HWP( 7 ) in Fig. [5) we identify 
(a) and (b) in the figure with the first and second spa- 
tial input modes). The two spatial modes are then com- 
bined at a polarizing beam-splitter (PBS in the figure) 
that reflects vertically polarized photons. The polariza- 
tion in one of the outgoing modes ((c) in the figure) is 
then flipped by a half- wave plate (HWP(7r/4)), and the 
spatial modes are split into two parts each, such that 
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there are four final modes which can carry two possible 
polarizations (in the figure: (e), (f), (g), (h)). The scat- 
tering matrix that describes the evolution of the four dis- 
tinct input modes (1-ff, IV, 2H, 2V) to the output modes 
(1H, IV,..., 4H,4V) reads 



W = 



V2 



/-s0-s0c0c0\ 
c c OsOsO 
1 1 
\ 00101/ 



(A10) 



where s = sin(27) and c = cos(27), i.e. the setup manipu- 
lates the polarization, and it depends on the polarization 
rotation angle 7. The result of a post-selected event with 
one particle in each output mode is a final state of the 
form 



fs - UV pulses 
Type II SPDC 




r^©^ ^©^L r±®~*L r*-®^ 



|<l 1<? ) = 



2V3 



sin^(2 7 ) |GHZ 4 ) 



'1,2 



(All) 

3,4) 1 



i.e. it is a coherent superposition of a bi-separablc double- 
Bell state and a four-particle GHZ-state, where 



IGHZ4 



1 ( 

2 

1 



\V,H) k ,), (A12) 



V2 



H, H, V, V) + \V, V, H, H)) .(A13) 



The success probability is (cos(47) 2 + 2 sin(27) 4 ) /12. 
The parameter 7 allows one to tune through several 
SLOCC-inequivalent entangled states of four photons 



Appendix B: Derivation of bounds on the 
dimensionality 



We derive here the upper bounds (37), (38) for the 



dimensionality dim(Sf) and dim(Se) of the manifold of 
accessible states, respectively. For this purpose, we re- 
group the matrix elements of W that appear in the state 
representation ( 30 ) in a smaller matrix W' of dimension 
N x dN: 



W' kd = W {k . 



■l)d+l,j> 



(Bl) 



such that W contains all relevant parameters (remember 
from Section II F 2 that we can always assume the initial 
separable state = |1, . . . , 1)). Formally speaking, in 
( 32 ) the dN 2 elements of the scattering matrix Wj k are 

mapped onto the ^-dimensional space of the gj 1 ,...j N ■ 
The dimensionality of the accessible manifold Sp is the 
effective number of complex parameters that are neces- 
sary to describe the 9j lt ...j N . 



FIG. 5: (color online) Creation of a family of entan- 
gled four-photon states. Two input modes (a) and (b) 
are combined at the polarizing beam-splitter (PBS), and they 
are subsequently split into four modes (e,f,g,h) by two addi- 
tional beam-splitters (BS). The initial state (A9i is thus scat- 
tered on a setup which implements ( |A10[ ) . The final quantum- 
information state |A12l depends on the polarization-rotation 
parameter 7. Courtesy of W. Wieczorek, C. Schmid, 
N. Kiesel, R. Pohlner, O. Giihne and H. Weinfurter |3T1. 



1. Fermionic qubits 

The dimensionality of the accessible space for fermions, 
dim(Hp), follows from mathematical properties of the de- 
terminant, since the coefficient 9j 1 ,....j N in (30) can be 
expressed as the determinant of sub-matrices of W' . 



We start with the case of qubits for simplicity, and 
re-write W as a vector of 2N, TV-dimensional column 
vectors, 



W = (wi,W 2 , ■ • -,W 2 n) ■ 



(B2) 



With ( 30 ) , we can write each coefficient of the state rep- 
resentation, §j lt ...,j N , as the determinant of the square 
matrix that consists of the respective columns of W , i.e. 



9h,...j N =Vet[(w jl ,w 2+j2 , . . .W2 N -2+ jN )} . (B3) 

If W does not have full rank (N), all coefficients (jj u ...j N 
vanish, because any set of N vectors w kl , ■ ■ ■ w kN is 
then linearly dependent. We can thus assume that W 
has full rank, and that the N vectors with odd in- 
dex, wi,W2+i, ■ ■ - W2{N-i)+ii are linearly independent. 
This assumption does not restrict generality, since sim- 
ple re-shuffling of the matrix and re-labeling of degrees 
of freedom can always provide us with such a situation. 
Consequently, the remaining vectors with even indices 
(1S2, W4,, . . . ) depend linearly on those with odd index, 
and we can always find coefficients c>^ k such that 



Vfc : w 2 k 



N 

£< 

3=1 



„(2(j-l) + l) 
"Ik 



W 2 (j-i) + i, 



(B4) 
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i.e. the vectors with even index can be expressed as lin- 
ear combination of the vectors with odd index. In or- 
der to express ( |B3[ ) in simpler terms, we use det(AB) = 
det(A)det(-B), or, equivalently, 



Det 




xDet[(u 1 ,u 2 ,...Ujv)]. (B5) 

In particular, the coefficient <7i,i,. ..1 is given by the de- 
terminant of the matrix consisting of the odd vectors 



wi , W3 , W5 . . . W2N- 1 (see ( 30 ) and ( B3 ) ) : 

<7i,i...i = Det [(wi,iM3, ■ ■ • w 2 n-i)} —■ V. (B6) 

Since the vectors W2k depend linearly on the W21+1, we 
can also express all other coefficients gj lt ...j N as a prod- 
uct of T> with determinants of matrices co mpo sed of the 



coefficients 
01,1.. .1 

02.1...1 

ffi,. ..,1,2 

32,2,1,...,! 



>0 



Therefore, i.e. exploiting (B5|, we have 



V. 



V ■ c. 



(i) 



V ■ Det 



(4") 



V-a 



(2JV-1) 
2N 



V ■ Det 



\^ C 2JV J 



2? • Det 



'2 c 2 
„(1) „(3) 



In particular, for all jk are odd, and 

■ ■ ■ Jn) is the empty 0x0 matrix: All rows and 
columns have been deleted, and g% t i ... t i = 2? - we set 
the determinant of the empty x 0-matrix to unity, for 
convenience. For 32, 2, 2. ..2, one needs to compute the de- 
terminant of the full matrix C = C(2, 2, ... ,2), since all 
jk are even. In general, we can rewrite the coefficients 
gj 1 ,...j N in the state e xpan sion as a product of two deter- 



minants, just like in (B5|: 

= V ■ Det [C(j 



■Jn)]- 



(B9) 



That is to say, the 2 N principle minors of C (includ- 
ing the empty matrix and the full matrix) are the co- 
efficients in the state representation. Since the values 
of wi, u?3, . . . w^n-i only affect the determinant T>, they 
contribute to only one free parameter. Consequently, at 
this stage, the maximal number of independent param- 
eters is bounded from above by iV 2 + 1 (iV 2 parameters 
for the matrix C and one parameter for T>). 

Given the functional form of gj 1 ,...,j N in terms of the 
principle minors of a matrix, we can actually find the ex- 
act number of independent parameters: The 2 N principle 
minors can be expressed by only N 2 —N+l parameters (a 
result which dates back to 1893 [92-494]). Consequently, 
we find 



dim(S F ) = N 2 -N + 2. 



(BIO) 



The principle minors also allow us to find a minimal 
explicit parametrization for accessible fermionic states. 
Since the set of all N (1 x l)-minors, all (N(N - l)/2) 
(2x2)-minors and all (N-l)(N-2)/2 (3x3)-minors that 
include the first row and the first column are sufficient to 
span the set of all 2^ — 1 minors [Mj , we can also express 
all accessible states in terms of these minors. 



2. Fermionic qudits 



3l,.. .1,2,2 — 



T) |J 2 ^-3)J2^- 
u I c 2N-2 C 2N 



_ (2JV-1) (2JV-3) 
c 2W-2 2N 



V ■ Det 



(27V-3) (2AT-1 

C 2JV-2 C 2N~2 

(2JV-3) (2JV-1 

C 2N C 2N 



(B7) 



For a general systematic representation of the 9j 1 ,...,j N 
as they are given in (B7), we group the coefficients 
„( 2 (j-i)+ 1 ) m a S q uare N x N matrix C with entries 



2k 



k,j 



(2(j-l) + l) 
-2k ' 



(B8) 



such that the state coefficient gj lt ...j N is proportional to 
some principle minor (or coaxial sub-determinant) of C 
[61] , i.e. the determinant of the sub-matrix C(j±, . . .j'jv) 
that is constructed as follows: For each k (with 1 < k < 
N), delete the fcth row and the kth column in the matrix 
C if jk is odd. The matrix C(j±, . . . jjv) is thus an m x m- 
submatrix of C, and m is the number of even coefficients 
jk, such that < m < N. 



In the case d > 2, we have a (d ■ N) x N dimen- 
sional matrix W' that determines the scattering pro- 
cess. The vectors w 2 , ■ ■ ■ ,w d ,w d+2 , ■ ■ ■ ,w 2d , ■ ■ ■ ,w Nd 
can be expressed as linear combination of the vectors 
wi,w d+ i, . . . ,uJ(at-i)<j+1) such that, in analogy to (B4|, 
we find 

N 



Vfc 7^ 1 + d ■ I : Wk — 



and the c^ 1 fill a matrix of size (d — l)N x N: 



E 



„(d(i-i)+i) 



Wd(j-1) + 1, 



(Bll) 



/ c 



C = 



(i) 
2 



„(d+l) (2d+l) 



„(1) (d+1) (2d+l) 



-d 
.,(1) 



(d+1) (2d+l 



„(<*+!) 



d 

„(2d+l) 
L d+2 



(1) (d+1) (2d+l) 
dN c dN c dN 



,((iV-l)d+l) 

,((i\r-i)d+i) 



,((iV-l)d+l) 

„((iV-l)d+l) 
"d+2 



J(N-l)d+l) 

-dN 



(B12) 
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In other words, for each additional dimension, we have 
to add N rows in the coefficient matrix. The d N coef- 
ficients gj 1 ,...j N correspond to d N minors of the matrix 
C. Consequently, the accretion of one dimension leads to 
up to TV 2 new parameters in the matrix C and thus to at 
most TV 2 new independent parameters. We thus have 



dim(H P ) < (d - 1)7V 2 -N + 2, 



(B13) 



which is the general bound for the dimensionality of the 
manifold of accessible states for TV fermions with d in- 
ternal states, when the fermions start in a non-entangled 
state. 



Bosons 



For bosons, the coefficients (30) are not the determi- 
nant of a sub-matrix of W , but they instead correspond 



to the permanent [68 70 . Since the permanent of a prod- 
uct of two matrices is not necessarily the product of the 
permanents of the two matrices, the argument that al- 
lowed us to reduce the number of parameters for fermions 
immediately breaks down for bosons. Consequently, the 
lack of algebraic symmetry properties of the permanent 
leads to a larger value for the upper bound to the dimen- 
sionality of accessible states for bosons than for fermions. 

In order to see this, we start with the dN 2 parameters 
W' k l that govern the d N coefficients gj 1 ,...,j N of the state. 
In the following, we identify 2N — 1 parameters that con- 
tribute to the normalization and to the global phase of 
the final state, and thus eliminate 27V — 2 of these. 

Firstly, a complex factor that is applied to each group 
of d columns leads to a global factor in the final state. 
More precisely, a factor ai to the Ith group of d columns 



corresponds to 



^j,(l-l)-d+k 



(B14) 



for 1 < j < N, 1 < I < N, 1 < k < d. By inserting the 
scaled W into (30), we see that these factors lead to a 



global factor that only affects the normalization and the 
phase of the final state, 



N 



9j u 



II"/ J 9h,...,j N - 
V/=l / 



(B15) 



In close analogy, the application of a factor to each row 
ofW, i.e 



also gives a global factor 



9jl,—,JN 




,3N- 



(B16) 



(B17) 



One of the N factors /L can be absorbed by the ctk, such 
that the a.j and (3k only provide 27V — 1 independent 
parameters. In contrast to the determinant, there are no 
further symmetry properties that allow us to eliminate 
other parameters, which leads to the following bound for 
bosons: 



dim(S B ) < dN 2 - 27V + 2, 



(B18) 



which is always equal or higher than the corresponding 
bound for fermions, Eq. (B13). 
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